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Distribution System Topology Detection Using 
Consumer Load and Line Flow Measurements 

Raffi Avo Sevlian and Ram Rajagopal 


Abstract —We present a topology detection method combining 
smart meter sensor information and sparse line measnrements. 
The problem is formulated as a spanning tree identification 
problem over a graph given partial nodal and edge power How 
information. In the deterministic case of known nodal power 
consumption and edge power flow we provide sensor placement 
criterion which guarantees correct identification of all spanning 
trees. We then present a detection method which is polynomial in 
complexity to the size of the graph. In the stochastic case where 
loads are given by forecasts derived from delayed smart meter 
data, we provide a combinatorial complexity MAP detector and 
a polynomial complexity approximate MAP detector which is 
shown to work near optimum in all numerical cases. 


I. Introduction 

The need for advanced controls in the distribution system 
is an emerging topic in power system and the controls com¬ 
munity. Proposed computational models for problems such 
as dispatching of distributed energy resources Q, El or 
coordinated voltage control 13,0, Q, la assume known 
system topology and network parameters. In reality customer 
level feeders are not known to within the accuracy as that of 
the transmission system. Therefore, situational awareness at 
the distribution system level is needed for future smart grid 
applications. 

The classical power systems state estimation approach has 
has been adopted by a number of authors. Such a method 
relies on the Generalized State Estimator (GES) Q, 0 which 
will solve for system voltages and breaker statuses. In a GES, 
breaker status variables are included into the state estimation 
equations as {0,1} variables that are usually relaxed and 
solved by rounding strategy El- These methods have been 
applied to radial distribution system with changing topology 
as well. There is recent interest around distribution systems 
topology detection. In ifTOl . a traditional GSE is employed in 
identifying the correct topology in a distribution system. They 
run a traditional WES state estimator and use dummy variables 
for breaker status indicators. 

Many non-conventional methods provide new insights into 
the estimation and detection problems. In Ea the authors 
introduce the use of high frequency micro-Phasor Measure¬ 
ment Unit (p,PMU) data in the topology detection task. They 
propose a method of comparing simulation and measured 


pPMU data for each topology. In ifT^ the authors develop a 
voltage time series approach to identifying topology changes 
relying on high accuracy /rPMU data. In ifTSll . lfT4l the 
authors develop a method using nodal voltage analysis to 
reconstruct the grid topology in real time. Similarly, na 
presents a method of reconstructing the network based on 
voltage correlation analysis. 

We should note some important features of each work. In 
m, im no approximations are made in measurements or 
system models and make no use of the tree structure of the 
distribution system. In ifTOl all forms of measurement are con¬ 
sidered since they rely on a general model formulation. In im, 
high frequency voltage magnitude and phase measurements are 
assumed at each bus. In El, El, M structural properties 
of the network are exploited, and linear models relating true 
power injections and observed voltage magnitude/phase are 
used. Additionally none of the method are used to provide a 
sensor placement strategy separate from heuristic evaluation of 
various scenarios. In El, the authors provide a method similar 
to El, where correlation analysis is used. However they 
rely on long time captures, so their method is more in line of 
network discovery not real time topology detection. In the 
authors present a general state estimator based method that is 
used in topology detection. They use the method for a heuristic 
placement strategy. 

Our work differs from these more traditional methods for a 
number of reasons. We focus on a simplified model measuring 
only power flow in a network. Our model does not consider 
voltage measurements or losses in the network. Second we 
make the assumptions similar to E3, EH where we aim 
to detect the operational tree structure of the network. We 
propose a general Spanning Tree Identification Problem and 
show the necessary and sufficient conditions for identifiability 
of all spanning trees under the measurement and power flow 
models. 

The paper is organized as follows. Section |I^ formulates 
the problem of topology detection. We present the load and 
sensor models used in the identification problem along with 
the sensor placement optimization. Sections ?? and|V]solve the 
detection and sensor placement problems in the deterministic 
and stochastic cases respectively. Einally, numerical demon¬ 


strations are given in Section VI 
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II. Problem Eormulation 

The topology detection problem is formulated as a spanning 
tree identification problem. Consider a power distribution 
network that has a planar graph structure but is operated as a 
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spanning tree at a given moment in time. Power is supplied 
from the substation to a set of feeders which are the at the root 
of multiple loads. In network reconhguration, sets of breakers 
and tie switches can reconhgure themselves such that all loads 
are connected and no feeders are connected. 

In the detection problem, we assume the breaker/tie switch¬ 
ing statuses are not available and must be inferred from the 
available power flow information. This is of practical concern 
since generally only feeder head breakers are automated, while 
other devices are mechanically set and reset by work crews. 
A usual concern in such situation is removing human error 
by automatically detecting the status of these devices from 
information currently deployed on the network. We assume 
the following information is available: 

1) Load pseudo-measurements in the form of forecasts 
which have forecast errors that must be taken into 
account. This information is available from Smart Me¬ 
tering Infrastructure. 

2) Real-time measurements of the power flows on a fraction 
of the lines obtained using noiseless sensors placed 
on the selected lines. These are modeled as noiseless 
since the errors of these measurements are negligible 
compared to those of the load pseudo-measurements. 

The detection and placement problem is solved in two 
scenarios: (1) Deterministic Case: loads and flows are known 
perfectly, (2) Stochastic Case: loads are known with uncer¬ 
tainty due to forecasting error. 

III. System Model 

A. Island Graph 


IEEE 123 Node Test Feeder 



together, while source of power are represented as feeder 
nodes. For example, sources 150, 251, 195 and 451 in Figure 
l(a)|are Feeder nodes Fi, F 2 , F 3 , F 4 in|l(b)| Additionally, the 


set of loads connected to these feeders, and each other are 
converted to loads vi^ 


, V 4 in Figure 1(b) 


In the island graph, we must maintain a switch configuration 
such that no feeders are connected to each other. The connec¬ 
tion of two feeders, for example by both wi and ws being 
closed will violate the tree structure of the network since two 
feeders will be energizing the same set of loads. 


B. Island graph 


TABLE I 

Mapping of vertices and edges to construct island graph. 


Island Graph 


IEEE Test Eeeder 

Switch 

Node 

Load ID 

Wi 

F2,Vs - Vi 

(250 - 251) 

W2 

Vi — V2 

(18-135) 

W3 

V2 - V3 

(151 - 300) 

W4 

Vi - ve 

(97 - 197) 

W5 

Vi - F3 ,Vt 

(450 - 451) 

We 

V3 - Vi 

(54 - 94) 

Wj 

Vi - V3 

(13-152) 

We 

Fi,ve - Vi 

(149 - 150) 

Wg 

Fi,ve - Fi,vg 

(95 - 195) 


To simplify the analysis of the switch configuration problem 
we modify the island graph. This is done by the following 
steps. (1) Loads from the island Graph are converted to vertices 
in the graph. (2) Feeders Fi.. .F 4 are turned into vertices. 
(3) Switches in Island Graph are converted into undirected 
edge in the graph. (4) A virtual root node and directed edges 
(vr —>■ Fj) for all feeders are added. 

A complete 


The result is the Island graph in Figure 2(a) 


mapping between the reduced feeder and the Island graph is 
given in Table For example, switch W 2 connects nodes vi 
and V 2 in the island Graph, or loads 151 and 300 in the test 
feeder. In the Island graph, it is now an edge between the 
vertices vi and V 2 - 


Fig. 1. 1(b) I Typical test feeder with sectionalizing switch operation. |(b)| 

Island graph simplifying topology of feeder. 

Consider, for example, the IEEE 123 node feeder in Eigure 
There are multiple loads connected to each other and 
separated by switches. We can then think of the network as a 
set of connected islands, separated by various switches. The 
typical feeder will have multiple power sources i.e. feeders. 
The set of switches are appropriately set so that each group 
of islands is connected to a load while ensuring no flow from 
one feeder to another. Therefore the distribution system must 
always maintain a tree structure. 

This network can be reduced to a simplified graph as shown 
in Eigure [T(b)] In this island graph, islanded loads are lumped 



Fig. 2. |(a)| Island graph showing potential topologies achievable. |(b)| Example 
spanning tree in network. 
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The Island graph G — {V, E) is the network used in the 
remaining analysis. Denote the added edges r = {e S i? : 
e = {vr, Fi) VFj}. This construction leads to a simple method 
for enumerating each valid topology of the Island graph. 

Consider T to be any spanning tree over G and T the set 
of all spanning trees constructed on G. We refer to the set of 
spanning trees containing subtree r as Tt-. Figure [^represents 
an example spanning tree that can be constructed. The number 
of spanning trees is indicated as |T| = T. 

Proposition 1: The set T,- represents all valid switch con¬ 
figurations. 

Proof: Any T € T,- will maintain that every island 
is connected. Additionally no feeders are connected besides 
through the virtual node. ■ 

Enumerating every spanning tree of the island graph which 
includes subtree r will generate the set of topologies corre¬ 
sponding to valid switch configurations. Note the naive method 
of enumerating all switch positions will enumerate all 
spanning forests of the graph which include the spanning trees 
associated with normal operation as well as the spanning tree 
under any number of outages where there are disconnected 
loads. In related work, the authors address the outage detection 
problem under the same measurement, detection and optimal¬ 
ity criterion in ifTTl . 


C. Load Model 

Each node excluding the virtual root and feeder nodes 
in the graph has a consumption load a;„. We assume for 
now that the load is single time invariant value. The pseudo¬ 
measurement of each load is Xn- We denote the pseudo¬ 
measurement error of load Xn by e„ = Xn — Xn- We assume 
that the loads are single phase real power quantities and 
the errors are mutually independent random variables: e„ ^ 
W(0, tJ^). Therefore, a;„ ~ A^(i„,cr^). x and x represent the 
vector of true loads and that of load pseudo-measurements; 
thus X ~ A^(x, S). Covariance matrix S being diagonal. 


D. Measurement Model 

Eor any edge e, we denote by s the power flow on it to 
all active downstream loads. The measured flow depends on 
unknown network topology and true loads. 

The sensor placement Ad C i?, is a subset of edges of the 
network. We assume that the magnitude and direction of power 
flow is measured. Additionally, we assume that the power 
flow measurements are noise free, with errors only with load 
pseudo-measurements. 

In general for a topology T € T and measurements 
s = {si,..., sm} the flow is Sk{x) = J2v,evUT,M) 
Where Vk{T,Ai) is the subset of downstream nodes for a 
particular placement and spanning tree. 

Given a general graph where all nodes can have some 
consumptions, we can represent the observations with 


s(r,x) = r(r,M)x 


( 1 ) 


where r{r,M) G {0,1}I'®I^I'^I. 

Eor instance the sensor placement in Eigure 2(a) (where all 
feeder nodes have 0 consumption) we have: 


1 ) sensor placement Ad = {ei, eg}. 

2 ) subsets Vi{T) = {vi,V 2 ,vz} and V' 2 (T) = {t^ 3 ,r' 4 } 
The observation matrix is then the following: 


r(r,M) 


110 0 1 
0 0 110 


( 2 ) 


E. Network Flow Representation 

In much of this work, we use the combinatorial graph and 
network flow representations to solve various the detection and 
placement problems. Relaxing the spanning tree assumption 
for now, we can represent network flows on a graph as follows. 



Fig. 3. |3(a)| Simple undirected graph. [3(b)] Network flows and direction (on 
edge), consumption and production of value (boxed) indicated. 


Given a graph incidence matrix B G {—1, 0, 
where each undirected edge has a pre specified direction: = 

(r’njt'm) on which to assign columns of B as follows. 

{ -fl if Ui is the originating node of edge 

— 1 if Ui is the terminal node of edge ej (3) 

0 else 

We then denote the flows on each edge of the graph as f G 
and the consumption vector y = [l^x — x]^. Therefore, 
the flow constraints can be represented as: 

Bt = y. (4) 


Eor example, the flow constraints for the undirected graph 
in Eigure |3(a) is the following: 


■ +1 

0 

0 

0 

0 

0 ■ 


1 1 

1_1 


yo 

-1 

-fl 

0 

+1 

+1 

0 



yi 

0 

-1 

-fl 

0 

0 

0 


= 

2/2 

0 

0 

-1 

0 

-1 

+1 



2/3 

0 

0 

0 

-1 

0 

-1 



. 2/4 . 


As an example, the following solution holds: f = 
[4 2 1 0.5 0.5 0.5]^ and y = [4 - 1 - 1 - 1 - 1]^ as 
shown in Eigure |3(b)| Note that negative fi indicate a flow 
opposite to the arbitrary direction specified in the incidence 
matrix. The sign of indicates a source (positive) or sink 
(negative) of energy. The graph indicates that vq is a single 
source for each sink in the network. 
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Assume a spanning tree T is encoded by the vector w, 
indicating the edges that are removed. Then the incidence 
matrix of a particular tree G will satisfy 

^wf = y- In the case of a spanning tree, the reduced matrix 
B^, where the first column is removed is an |y| — lx |T^| — 1 
invertible matrix. We can represent the matrix r(T, Al) by 
the following procedure. 

r(r,M)=A(M)B:;-^ ( 6 ) 

Matrix A(A4) indicates the edges that are being measured, so 
A{k,mk) = IjVmfe G Ai. The superscript i?’’ indicates that 
the hrst row the matrix has been removed. 


IV. Deterministic Case 

A. Detector Operation 

Given nodal consumptions x and observed flow s the 
following deterministic program is proposed: 


flndf,w (OPT-1) 

s.t. 

\fi\ < |x|iu* (7) 

B{ = y (8) 

fi = s, Ve, G EiM) (9) 

w G SP-TREE (10) 


The feasibility program is mixed integer in with boolean 
vector w G {0,1}I^I for the edges in G as well as the flow 
along each edge f G The network flow constraint is 

given in 0. Constraint 0, limits the edge flows to either be 
set to zero, or be fully unconstrained. Constraint 0, sets each 
observed edge to the sensor value while ( [TOl i constrains the 
status of edges in G to form a spanning tree. 

r = {rGT:s = s(r,x)} (ii) 


The naive detector in o will enumerate every spanning tree, 
then given the nodal consumptions, evaluate the theoretical 
flow value s(T) and compare it to the observed flow. There 
are two questions regarding this simple technique: 


1) Correctness: How to guarantee that this method will 
return a unique and correct spanning tree. How to 
guarantee that each spanning tree results in a unique 
flow observation. 

2) Efficiency: Enumerating each spanning tree then per¬ 
forming a comparison will lead to a detector of 0(|T|) 
complexity. Typically |T| grows exponentially in the size 
of the graph. 


The question of correctness is solved in IV-B and an 
efficient polynomial complexity algorithm is introduced in 

USD] 


B. Detector Objective: Spanning Tree Identifiability 

We only require a sensor placement such that different 
spanning trees lead to unique measured flow values. The 
following definition is of use: 

Definition 1: The set T is identifiable if VT, T' GT where 
T f T' v/e have that s(T, x) f s(T',x). 


In the deterministic case, we desire a placement M. such 
that T is identifiable. A naive method of determining whether 
T is identifiable is to enumerate every single spanning tree, 
generate r(T, M) to verify that for no two trees T{T, Ai) = 
^(T^Ad). This naive procedure has 0(|Tp) complexity, 
which can be quite intensive since |T| is generally very large. 

C. Deterministic Sensor Placement 

An important result is the following: 

Theorem 1: T is identifiable if and only if the graph E{G)\ 
E(Ai) of the Island graph forms a spanning tree. 

We refer to any placement Ai satisfying Theorem [T] as a 
valid placement. 

The proof relies on graph theoretic arguments as opposed 
to linear algebraic arguments on the relation in eq. 0. The 
intuition of Theorem is that to have observability of all 
spanning trees, we must have a sensor placement such that any 
cycle that can be constructed on G will have some flow sensor 
on it. This corresponds to the dimension of the cycle space, 
referred to as the circuit rank p,{G) = \E\ — \V\ + I, which 
is the minimum number of measurements needed to correctly 
detect all spanning trees on G. This gives us a 0{E) verifiable 
condition to ensure that all spanning trees are identifiable as 
opposed to C)(|Tp) with the naive method since we only need 
to remove particular edges of the graph and test if it is a 
spanning tree. 

Theorem provides a way to construct the set of all 
placements where identifiability is achieved. 

Eirst consider the function h{T) = E{G) \ E{T), which 
merely extracts the edges in G that are not in T. An obvious 
consequence to Theorem is the following: 

Corollary 1: The function h : T —> AI is a bijection 
between the set T and the set of all valid placements M. 

The following is useful: 

Remark 1: Corollary implies that |M| = |T|. 

Corollary [T] is quite important from a placement perspective 
since it actually yields a method to generate a valid placement 
in the deterministic case. Also, it allows us to enumerate all 
valid placements for a graph. This is important when dealing 
with a stochastic case where sensor placement relies on mostly 
evaluating each placement in M. 

In the case of valid placements on the Island graph this 
set is restricted, since having sensors on edges in r, since 
they have no physical meaning. The restricted set is given by 

Mr = {fiT)\T GTr}. 

D. Spanning Tree Detection via Relaxed Flow Solution 

Theorem [^states that if E(G)\E{M) forms a spanning tree, 
there is a unique observable flow for each unique spanning 
tree. However, it provides no efficient method of actually 
decoding the correct spanning tree. 

f* = {f : St. B f = y and fi = Si Ve^ G E{M)}. (12) 

We show that solving a combinatorial search over all spanning 
trees, reduces to solving a relaxed form of the problem, over 
the fully connected graph. Surprisingly, the solution to the 
linear equation over f G recovers the sparsity pattern 
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in f corresponding to a spanning tree without any sparsity 
inducing heuristics. 

We can represent the conservation of flow in eq. 0, by 
partitioning the incidence matrix B and flow vector / into 
observed {Bm 1 ^m) and non-observed j fjv) components. 

Therefore, we can represent the network flow with the 
following; 

Dr Dr 

0 I 

Lemma 1: For the sensor placement condition in Theorem 
[^where E{G)\E{M) is a spanning tree, rank(Bjv) = N—1. 

The matrix Bm G is of rank N—1. The sub 

matrix B"^ with the first row removed is invertible, which 
leads to the following solution to ( [T2l i. 

f(x,s) = (B]v)"^(x-B];^s) (14) 

Where we have removed the first row of Bm and Bm- Next 
we must show that the solution to this is in fact the correct 
spanning tree on the graph. 

Theorem 2: If Ni satisfies the condition in Theorem [T] the 
solution vector f = [f(x,s) s]^ encodes spanning tree T- 
The search over the set of spanning trees, can be replaced by 
solving a set of linear equations eq ([T^. 




X 

fM 


s 


(13) 


V. Stochastic Case: Detection 
A. Maximum a Posteriori Spanning Tree Detector 

Recall that in the stochastic case we integrate (1) error 
free flow information on edges (2) load pseudo measurements 
at nodes. Combining these two leads to a hypothesis testing 
problem as follows: Given the forecasted loads x, the true 
loads at each node are given as: x ~ 7V(x, cr^J). Therefore, 
under a particular hypothesized spanning tree T, the true flow 
would be distributed as: 


s \{%, x} ^ iv(r,x, a^r.rf) (i5) 

(16) 

= s(7^,x) + (17) 


Each line is an equivalent definition. The matrix F^ in eq. 15 


denotes F (T, M) for tree % under a known, fixed placement 
The term s(7i,x) in eq. 16 17 indicates the theoretical 


flow that should be observed under the forecast of nodal 


consumption. The error Cs.i in eq. 17 is a zero mean multi¬ 
variate gaussian with covariance matrix Assuming uniform 
prior on hypotheses, the general MAP detector reduces to a 
maximum likelihood detection of each spanning tree: 


T = argmaxPr(s | x, %) (18) 

Ti GT 

The MAP detector is optimal in minimizing the mean square 
error but suffers from the 0(|T|) complexity of the determin¬ 
istic case. 


B. Flow Based Approximate MAP Detector 

The flow solution in eq. can be used to construct 
an efficient hypothesis detector which does not have 0(|T|) 


complexity. Recall in the deterministic case, the flow is f (x, s). 

In the stochastic case, s and x is given, therefore computing 
f(x,s), can be useful. The goal is to compute the flow given 
the measurements available (s and x), and use it for detection 
of the correct spanning tree. First, we can determine the 
distribution of this flow vector conditioning of a candidate 
hypothesis %\ 

f(x, s) I {71} = x)) (19) 

= + ( 20 ) 

= ( 21 ) 

-- W(f(x,s(7},x)),E/,j) : Ti is true (22) 

The LHS of eq. ([T^ is the distribution of the ’noisy-flow’ 
conditioning on a particular hypothesis 7). The RHS evaluates 
the flow equation ([T^ using the forecasted consumption x 
instead of the true value x. Finally, in (|2T]l-(|22ll the relationship 
between s(7},x) and s(7i,x), leads to eq. @-([n]i. 

The vector f(x, s(7},x)) is the flow from spanning tree 7) 
and nodal injections x. The true ’noisy-flow’ is distributed 
around this value, depending on the hypothesis T- 

A possible hypothesis test is the following; 

T =argmax Pr(f(x,s) | f(x,s(7},x))). (23) 

VTiGT 

This is no better than ( [TSl l for the following reasons: 

1) We must compute the hypothesis mean f(x,s(7i,x)) 
under every spanning tree %, it is still of 0(|T|) 
complexity. 

2) The covariance matrix is of rank p{G) (the rank of 
Fi) and not E (the size of f). Therefore EtJ is not 
positive definite and therefore the inverse { cannot 
be computed. 

An alternative is to test p(G) elements of the ’noisy-flow’ 
f(x,s) under the hypothesis that their true value is zero 
corresponding to the zero’s of the hypothesized spanning tree. 

Using the following shorthand: T = f(x, s) and fj- — 
f(x,s(T, x)) we can represent the variables in ( |2^ . Under 
hypothesis T, fr(0 = 0 for some set of indices i G Ij- 
where Ij- = (j : G E{G) \ E{T)}- 

Therefore, under a hypothesis T, we can calculate the 
following likelihood; 

f = argmaxPr(fo(ii),...,fo(i^)|fr(ii) = 0,...,fr(v) = 0). 

It-TgT 

(24) 

In this case, the reduced covariance matrix is potentially 
invertible. The following theorem relates the combinatorial test 
to a test of zero flows on the empirical flow. 

Theorem 3: The zero flow hypothesis detection in ( |24l i and 
the combinatorial flow hypothesis test in eq. are equiva¬ 
lent for A4 satisfying the placement condition in Theorem [T] 
of size p.{G). 

Therefore the MAP detector over the set of all spanning 
trees reduces to testing whether an observation vector Unfor¬ 
tunately we still need to enumerate |T| hypotheses. 

Fortunately, however vector form allows us to very effi¬ 
ciently prune out all but a few alternative hypothesis to test. 
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Intuitively, f(x, s) will have some very few small values which 
actually encode potential spanning trees, and many very large 
values which can just be pruned. 

This leads to the following two step approximate MAP 
detector: 

1) Evaluate the empirical flow f(x,s) via ( [T^ . 

2) Compute the minimum weight spanning tree solution on 
graph G where edges are weighted with — |f(x, s)|. 

C. Cycle Descent Approximate MAP Detection 

An alternative approximate MAP detector is based on 
generating single cycle edge exchanges Aci = {ci e'}. 
For ever edge in in the co-tree of current tree, the perturbed 
edge is chosen to be an edge along the fundamental cycle of 
Ci. The complete algorithm is given as the following: The 


Input: [1] Observed flows s. 

[2] Load Forecast x and Error Covariance E 

[3] Graph G 

Output: MAP Detection Hypothesis T 

1 // Must start with feasible Tree 

2 T ^ feasible-tree(so, x, G) 

3 FG ^ generate — fundamental — cycle (T) 

4 while Aloglik 7 ^ 0 do 

5 for Cfc e FG do 

6 loglik, T ^ local — update(T, Cfc, s, x) 

7 FG update — cycles (T, FG) 

8 end 

9 end 


cycle descent algorithm performs as follows: 

1) feasible-tree .• A feasible tree is generated which 
satisfies the pattern of observed zeros as the observed flows s. 
The edges of G are given weights as follows: 

• Edges with measurements are weighted as follows: zero 
measurements have weight 0. 

• Non-zero measurements are weighted with some positive 
value K. 

• Remaining edges are assigned a very large weight (« 
\E\ K). 

Finally, a maximum weight spanning tree is calculated on the 
weighted graph. This ensures that the initial spanning tree To 
satisfies I{s > 0} = I{s(T, x) > 0}. The remaining spanning 
tree is feasible, where evaluating Pr (s | Ti, x) will lead to a 
finite (albeit small value). 

One would assume that the maximum spanning tree output 
of FMST method can be used as a suitable initial tree. 
Unfortunately, the output of the FMST is not guaranteed to 
produce a tree that satisfies the observed non-zero flow. Edges 
that measured zero flow will not be included in the initial 
condition tree, however we are not guaranteed to include all 
edges with non-zero measured flow. 

2) local-update .• For cycle c^ and edge Ck, the next edge 
is chosen by evaluating the likelihood of function T—Ck+ej, 
where ej G E{ck)- Not every The next edge Cj is chosen 
which maximizes the likelihood. 


3} update-cycles .• After each edge exchange operation, 
the cycles must be updated. A queue is maintained for the 
edge Cfe to be processed, where the elements are updated. 

D. Local MAP Search 

A simple local search procedure is the following: 

1) An approximate detection is performed: (using FMST 
or Cycle Descent) with output T*. 

2) A subset of spanning trees ’close’ to T* are enumerated 
with likelihood function evaluated. 

We define trees near T*, to be all spanning trees that map to 
the same fundamental cycle basis. Enumerating the candidate 
edges of this set is of complexity 0{\E\p,{G)) and is discussed 
in in. 

VI. Numerical Experiment 
A. Deterministic Placement 

We test the placement problem on a set of planar graphs, 
shown in Figure]^ In both graphs, designate a single vertex as 
the source, which is indicated in the dashed horizontal lines. 
Graph Gi has Uroot = V 4 and G 2 has Uroot = 'Pi- 



Eig. 4. Sample graphs used to benchmark deterministic and stochastic cases, 
Gi with Uroot = V 4 ((Eigure [4(a)l . These graphs have high circuit rank but 
smaller |E|//i(G), therefore difficult to visualize stochastic placement. Graph 
G 2 (Figure |4(b)j used to illustrate stochastic placement as well as provide 
sub modularity counterexample. 


TABLE II 

Deterministic Topology Detection. 



KG) 

|T| 

mean 

e' 

std. 

max/min 

\E\/KG) 

Gi 

5 

391 

56.9 

26.7 

299/6 

2.8 

G 2 

5 

830 

139.3 

72.2 

185/10 

3.2 


To test the placement problem, we enumerate the set of 
spanning trees for each of the graphs. The method relies on the 
backtracking method developed in lfT9l . The simulation was 
implemented in MATLAB and deemed correct by checking 
that each spanning tree was unique and the number of test 
trees corresponded to those calculated from the matrix-tree 
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theorem EOl . where T = det(L„). Where Ly is the v minor 
of the laplacian matrix with the result being invariant to v. 

For the three graphs in Figure the graph statistics and 
experiment results are shown in Table [I^ We evaluate the 
experimental error rate 

£ = ^ (25) 

From Theorem [T] the missed detection error must be zero. The 
computed e was zero in both cases, as was expected. 

We evaluate the output according to an arbitrary input 
because we would like to compare it to the case where only 
magnitude and not direction is measured. This is a common 
type of power system measurement, where no phase is reported 
and power direction flow is unknown. In this case, we evaluate 
e' which now compares |s(7i,x)| ^ |s(7),x)| instead. The 
computed values for e' are shown in Table [1^ We evaluate 
each valid placement in M to illustrate the importance of flow 
direction. The value reported in Table is the mean missed 
detection ± the standard deviation. This verihes that different 
placements result in different unsigned missed detection rates. 
We see that if the direction of flow is not known, around 10% 
of the spanning trees are indistinguishable on average. 


B. MAP Detection Performance 



(a) (b) 


would make one naively guess that the missed detection error 
should not depend on any single placement. Within a graph, 
it is observed that the the mean length of all fundamental 
cycles associated with a placement is slightly correlated with 
the mean error. 

C. Approximate MAP Detector 



^ -5 


6i-io 

o 


o 

E 0 


>-1 


x10 

Greedy Failure 


I7"“ “ 

A" 


2 

4 6 8 10 

xio" 

Greedy Success 



5 10 

Iteration 


15 


(c) 



(d) 


Fig. 6. Comparison of mean missed detection error for optimal MAP 
detector, flow based approximate detector (FMST/FMST-MAP) and cycle 
descent detector (GR/GR-MAP) for the graphs G1 j5(a)l , G2 j5(b)l with 
respect to a. Example of cycle descent success and failure for two 

spanning trees. W)) Evaluated missed detection probabilities for the set of 
placements on the network Mr- For half of the placements, the conditional 
pdfs are far enough that there is negligible error. 


Fig. 5. Mean missed detection error for the graphs G1 j5(a)^ , G2 j5(b)| with 
respect to cr. 


The performance of the MAP detector is evaluated for each 
of the Graphs in Figure The one shot detector performance 
is evaluated with a uniform load mean of /i^ = 1 and forecast 
error of a. Figure shows the mean missed detection error 
over all hypotheses with respect to cr. 

A number of important observations can be seen from 
this analysis. Different graphs experience widely different 
behavior. For example Gi has many very short cycles where 
many of the spanning trees will see multiple zeros while G 2 
has only two cycles with high edge count per cycle. The mean 
errors in Gi are must smaller (roughly half) as that of G 2 . 

The sensor placement in the stochastic case has a dramatic 
impact on the mean missed detection error. This is slightly 
counter-intuiive, since a single placement which maps to a 
single spanning tree must correctly decode all spanning trees 
with low error. A symmetry between placement and tree’s 


The two approximate MAP algorithms are tested on Gi 
and G 2 , where the performance are shown in Figure 6 (a)|6(b) 
respectively. 

1) Flow-Maximum Weighted Spanning Tree (FMST): In the 
near noiseless case, the flow based approximate MAP detector 
performs identically to the combinatorial map detector, in 
both Gi and G 2 . This is because the ’noisy-flow’ values are 
very close to their correct values of zero. For this reason, 
the maximum weighted spanning tree algorithm can correctly 
reconstruct the correct spanning tree. 

In the high noise case, the algorithm fails worst in the 
high cycle count graph Gi, where the maximum spanning 
tree graphs very rarely match with the maximum likelihood 
output. In G 2 , however the two are nearly identical since the 
tree does not involve many Subsequent local search leads to 
an improvement in the mean missed detection rate. 

2 ) Cycle Descent Algorithm: The cycle descent algorithm 
has a very different performance than the FSMT algorithm. 
For both graphs, the performance is similar, in that a small 




































subset of trees in Gi will always fail regardless of SNR. For 
the remaining spanning trees, the algorithm converges to the 
MAP output. In the case of Gi, only 4% of the trees lead to 
a failure of the algorithm, for the used placement. Different 
sensor placements lead to different values of non-convergence. 


D. 123 Test Feeder 


As presented in Section we introduced a method of 
transforming the IEEE 123 test feeder into a reduced feeder 
(Eigure l(b)| l and an island graph (Eigure 2(a)| l. The results 
we present allow us some insight into the topology detection 
problem. Erom Theorem [T] the minimum number of flow 
measurements is /J.(G') = 4. Since the actual set of allowable 
trees is T,-, this restricts the set of trees to search over, in 
evaluating for each sensor placement: 

1) Mean missed detection error over all possible spanning 
trees: 


gi{M) = Y,Vr{T)Vr{f^T\T-M)- (26) 

Tgt 

2) Maximum missed detection error over all possible span¬ 
ning trees: 

g 2 (Af) =m|xPr(r ^r|r;A4). (27) 

Eor simplicity assume a single phase flow of only real 
power. Given the three phase complex power consumption of 
each load, we reduce this to a single value power value. This is 
the case, if the measurement on a line sum and report the sum 
of the three phases. Eor a rough estimate of the forecast error, 
we follow the forecasting model in ED, where day ahead 
forecasts are computed for varying aggregation levels. The 
forecast errors are used to construct the following scaling law 
for coefficient of variation: 

- = y^+41.9. (28) 

fi \ W 

. Since the loads of each individual island is quite large and 
beyond the critical load reported in ED, the CV of each island 
is close to the irreducible error of 6.3%. Eigure [6(d)] shows the 
performance of an initial stochastic placement for the IEEE123 
test feeder. We evaluate the set of restricted placements M^- 
and spanning trees Tt-, where |Mt-| = 44. Notice that for 
almost half of the placement the maximum error is negligibly 
small. One way to consider this is that although the various 
hypotheses means have signification error variances, some 
placements separate out the hypotheses sufficient such that 
they are all distinguishable. 

This example allows us to comment on the common under¬ 
standing of generalized state estimation. In current practice, 
determining the boolean values of breakers or In practice the 
set of all topologies generally map to range of ±5%. There¬ 
fore, when factoring in uncertainties, the missed detection 
rates can be high. On the other hand, measuring flows lead 
to very large changes in the observation vector over the range 
of hypothesis. 


VII. Conclusion 

We formulate the problem of detecting switch configurations 
on residential feeders as a spanning tree detection problem on 
an ’Island graph’. The detection problem relies on power flow 
measurements on edges as well as load information at nodes. 
The deterministic case leads to notion of tree identifiability 
and sensor placement conditions to ensure identifiability. Eor 
the stochastic case, we propose a greedy algorithm which leads 
to a locally optimal solution. 
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Appendix 


G = {V,E) Undirected Graph G; Vertices U; Edges £. 

r 

Spanning tree on G. 

T 

Set of all spanning trees that are constructed 
on G. 

T 

Set of edges in constructing Island graph. 

M 

Sensor placement {A4 C £). 

M 

Set of all sensor placements leading to identi- 
fiably in T. 


True and forecasted load of node 


Forecast error variance and covariance matrix. 

nr.M) 

Observation matrix for tree T, sensor place¬ 
ment Ai. 

^obs 

Set of measured power flow. 

s(T), s(T) True and predicted flow measured under hy¬ 
pothesis T- 

Uj 

Acceptance region of 7} over 7}. 

Ri 

Acceptance region of hypothesis hi over all 
alternatives. 

gi,2{M) 

Shorthand notation for (1) maximum / (2) mean 
missed detection probability over all hypothe¬ 
ses. 

c 

Cycle in graph G. 

C{G) 

Cycle space of G or set of all possible cycles. 

FC{T) 

Fundamental Cycle Basis of G constructed by 

r. 

FCm 

Fundamental Cycle Basis constructed by sen¬ 
sor placement AA. 

Afc 

Afc is A:* cycle in FCm- 

p{G) 

Circuit Rank of Graph. 

n{G) 

Number of connected components. 

AE 

Edge exchange operation to generate new tree: 

T^V- 

IC{c) 

Cycle-Sensor Map indicating all sensors on 
cycle c. 


A. Useful Graph Theory Definitions and Results 


Vector addition is defined as c' = ci © C 2 where new cycles 
are constructed via symmetric difference operation on edges: 
Ei®E 2 = {El U £'2) \ {El n £2)- The cycle space C{G) 
of the graph is the vector space of all possible cycles in a 
particular graph. 

For example consider the cycles in the graph G in Figure 
|7(a)| (assume dashed lines now solid): 

Cl = {67, Cl, 62, eg} 

C 2 = {eg, e 2 , eg, eg, eig, 612 } 

Cg = Cl © C 2 = (ey, ei,eg, e 4 ,eig,ei 2 } 


It is easy to see that ci, C2 and cg are all cycles in G. 

Circuit Rank: The circuit rank of a graph is given by p,{G) = 
|£| — \V\ + n{G) where n{G) is the number of connected 
components of the graph. For example the Island graph in 
Figure 2(a) has n{G) = 1, and /r(G) = |13| — |10| + 1 = 4. 

Cycle Basis: The analog of a vector basis for cycle spaces is 
the cycle basis. A basis Be C Be is the smallest number 
of cycles whereby all other cycles can be constructed via 
symmetric difference operations. The dimension of C{G) and 
Be is b{G), the circuit rank of the graph. Therefore, p(G) is 
the smallest number of cycles required to produce all other 
cycles on a graph. We use ^(G) and p, interchangeably, 
whenever convenient. 

For example, all cycles in the graph in Figure 2(a) can be 
constructed from p{G) independent cycles which form the 
cycle basis of the graph. 


Fundamental Cycle Basis: A Fundamental Cycle Basis FC is 
a cycle basis constructed using the following procedure: Given 
a spanning tree T, enumerate the set of edges in G but not 
T. Then for each e € £(G) \ £(T), construct T + e then find 
the single cycle c associated with e. So we can generate p{G) 
cycles in FC, which is the dimensionality of the basis. An 
equivalent definition for a Fundamental Cycle Basis is that 
each cycle will have one unique edge which is in no other 
cycle. 



(a) (b) 


Fig. 7. Two spanning trees F j7(a)l and F' l |7(b)k with sensor placement 
M. 

Refer to ll20l for a more thorough presentation. 

Cycle: A cycle c = {ei,...,eAr} is a connected subgraph 
where each vertex is of degree 2 . 

Cycle Space: The set Be is, the power set over the edge set 
Be = (0, Any cycle c is a vector in the space Be- 


TABLE III 

Fundamental Cycle Basis erom T, T' and cycle-measurement 
MAP /C(c) 


FC{T) 

JC{c) 

Cl = (ce, 62 ,Cl, 67 } 

/C(ci) = {si} 

C2 = {ei2, ei 3 , es, 64, ci, 67} 

A{C2) = {si,S 2 } 

C3 = (cs, 63,64, es} 

E-{cfj = { 52 , 54 } 

C4 = {69,68,64,61,67,611,610} 

^(^4) {'^11 *^2 ; '^4} 

FC{T') 

/C(c) 

Cl = {67, 66, 62,61} 

JC{ci) = {51} 

62 = {613, 68, 64, 62,66, 612} 

E,{c 2 ) = { 51 , 52 } 

63 = {65, 68, 64, 63} 

E-{cf) = {52,54} 

64 = {69,68,64,62,66,611,610} 

^“(64) = {51, 52 ,54} 


From spanning trees T and T' in Figure 7(b) 7(b) 


we 


construct the following Fundamental Cycle Basis in Table III 


Note that FC{T) f FC{T')- However, this does not always 
occur, see HH for more details. 
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B. Proof of Theorem I 


Recall, that theorem states that as long as E{G) \ E{M) 
forms a spanning tree, then for any two trees T, T', s f s'. 

We now prove Theorem [T] in the following steps. 

1 Introduce a cycle-measurement mapping object /C(c) 
which tracks sensors on a cycle and show that an in¬ 
dependence property if E{G) \ E{M) forms a spanning 
tree. 

2 We construct an edge exchange procedure which encodes 
the transition; T ^ T between any two spanning trees 
as a set of single cycle edge exchanges on the cycles of 
EC{T). We show that this encoding always exists. 

3 We show that sensor measurements in JC{c) decouple 
from one cycle to another under single edge exchanges. 

4 We use the independence of JC{c) and decoupling of sin¬ 
gle edge exchange measurements in c to show inductively 
that no multiple edge exchanges of any size exist where 
As = 0, if E{G) \ E{M) is a spanning tree. 


1) Cycle-Measurement Map: The cycle-measurement map 
encodes which measurements lie on the edges of a particular 
cycle. The map /C(c) is dehned according to cycles c G EC, 
for any arbitrary EC in G. 

Definition 2: With respect to some EC, a cycle- 
measurement map is /C : c —^ for all c S EC 

where Sk G A(c) if Sk is on an edge in c. 

We can also write it as A(c) = cfl E{Ei), though this is an 
abuse of notation. 

For the fundamental cycles associated with the tree in Figure 


7(a) 7(b) we construct the following map shown in Table III 


We now aim to develop some useful properties of this mapping 
function. 

Consider placement f4, and constructed tree T = E{G) \ 
E(JCi). We must have JC{ck) = {s/c} where Ck is the fc* cycle 
in EC{E). This is obvious by construction. 

We denote Xk to be the fc* cycle in EC{E{G) \ E{M.)) 
{ECm shorthand). From before, IC{Xk) = {sfc}- 

We see that by looking at /C(c), for any cycle in an 
arbitrary EC, we can encode the cycles construction using 
basis ECm ={Ai, ..., A^}. 

Lemma 2: If E{G)\E{Ei) forms a spanning tree, then for 
any EC and c G EC, c= 0 A^. 

fe:sfeG/C(c) 

Proof: Any cycle c can be represented as a combination of 
cycles A G ECm since the ECm is a cycle basis. If any other 
A' outside of the set {Afc : Sk G A(c)} is used to construct c 
then JC{c) will include edge containing a s'. Conversely if any 
additional X' is needed to construct c, it’s s must be in /C(c). 


Proof: Suppose that there exists some EC and partitions 
A, B where the terms ICa and Ais are equal. Since 1C a and 
ICb encode some cycle we have that ca = Cb- However, since 
CA and Cb are dehnition fundamental cycle’s in EC there will 
exist an edge in ca that is not in cb thus ca cb- ■ 

To see an example of this consider the cycle-measurement 
map generated by EC{E) in Figure 7(a) in Table III Notice 
that we cannot construct any partitions A, B where all the 
covered measurements are equal. For example if A = {ci, C 2 } 

and B = { 03 , 04 }, we have that IJ lC{ck) = {si,S 2 } and 

keA 

U = { 51 , 53 , 54 }. 

k£B 

Notice that the cycle-measurement-map 1C constructed for 
the network in Figure [2(^ does not satisfy that E{G)\E{M) 
is a spanning tree. In this case we have that ICicf) = {si}, 
lC{c 2 ) = {si, S 2 }, lC{cf} = {S 2 } and ^( 04 ) = (si, S 2 } for the 
same cycles. The cycle partition A = {oi}, B = {04} clearly 
leads to the independence property not holding. 

This result leads to the following equivalent results which 
are used in Section lB4l 

Corollary 2: A special case is that Vo, A C EC, 3s G /C(c) 

s.t. s ^ IJ /C(c), for any A not including o. 

k^A 

Remark 2: The subset G C EC, where \G\ = N will have 
at least N unique sensors. 

Remark 3: The indicator vector associated with each /C(c) 
are a set of linearly independent vectors. 

2) Edge Exchange Operator: We encode the transition 
between any any two spanning trees and show that any 
such transition can be represented as a set of single edge 
changes. To motivate this, consider the trees T and T' in 
Figure |7(a)| |7(b) The removed edges from G in each tree are 
Er = { 66 , 612 , 69 , 65 } and E'^ = { 67 , 613 , 69 , 65 }. In both 
trees, 69 and 65 do not change. 

The main question we want to answer is how to encode the 
transition between trees by single edge exchanges. Namely, if 
we define a AE operation, do we encode as AE = {cq ^ 
67, 612 613} or AE = {66 -)► 613, 612 -)► 67}. 


TABLE IV 

Fundamental Cycle Basis erom T, T' and cycle-measurement 
MAP /C(c) 


EC{T) 

cDEr 

cCEr 

Cl = {ee, 62, 61,67} 

{ee} 

{er} 

62 = {ei2, 613, 63,64, 61,67} 

{ei2} 

[er, eis} 

63 = {65,63, 64,63} 

{eg} 

{eg} 

64 = {69, 63, 64,61,67, 611,610} 

{es} 

{ee} 


Now we can prove a general case of ’independence’ between 
any two subsets of cycles and the measurements that are placed 
on them. 

Lemma 3: If E{G)\E{A4) forms a spanning tree, then for 
any EC and subsets A f B of cycles in EC we must have 
that: 

[JlC{ck)f\JlC{ck). (29) 


k^A k^B 



This can be resolved if we look at edge exchanges with 
respect to EC{E), as shown in Table IV Column 1 repeats the 
cycles in EC{T). Column 2 maps edges in Er onto EC{E) 
and Column 3 maps edges in E'j^ onto EC{E). This can be 


seen as an identical mapping function as IC{c) in Section B1 
except we replace edges with measurements on them, with 
edges in Er. Notice that we can now define a cycle by cycle 
set of edge exchanges, that dehne the transition from T —> T^ 
So on cycle ci we have Ae = {cq -a 67) and on cycle 62 we 
have Ae = (612 —>■ 613). 



















II 


Definition 3: An edge exchange with respect to J^C is 
l^E = {Aei,..., Ae^} where Ae/c = (cfc e'^.), e/c S Eh, 
e'f. e E'^ and ek,e'^, G c^. 

We can generate an edge exchange encoding as follows. 
First assign the edge in Eh used to construct Ck G EC{T). 
For e'j, use the following procedure: 

1) Starting at ci, set to be any element in cnE'^. 

2) For the Ck, set to be any element in Ckf^E'^, that 
has not already been assigned to previous e'^..., e'k-i- 

Lemma 4: For any two T, T', at least one edge exchange 
procedure always exists. 

Proof: The set of edges cr\EH are equivalent to the map¬ 
ping K{c). Therefore, Corollary]^ holds for each incremental 
cycle to be processed: that is, every additional c(^Eh set will 
have a one new edge that can be assigned to e'^. ■ 

Therefore between any two spanning trees there is a well 
defined set of single edge exchanges performed on the cycles 
of FC(T) which encode any arbitrary T ^ T' ■ 



Fig. 8. |8(a)| Illustrative example to see properties in Lemma 


4) Inductive Proof: Note our effort is to show sufficiency, 
for necessity we only need to consider the cycle c in Figure 


8 (a) If no sensor exists on q, E[G)\E{Ai) is not a spanning 


tree and every edge exchange leads on Ci will lead As = 0. 

From Lemma we can encode any T ^ T' transition as a 
set of cycle based edge exchanges. We now show inductively 
that no edge exchange exists which leads to As = 0. So T —>■ 
T' always leads to As 0. 

We show that given any number of non-trivial edge ex¬ 
changes where at least one Si s.t. As^ 7 ^ 0 . 

Base Case Assume c C EC, contain a non-trivial edge 
exchange. From Lemma |^(P2), Vs G JC, As 7 ^ 0. 
Inductive Hypothesis Assume multiple cycles C C EC, 
contain a non-trivial edge exchange, where \C\ = N. 
Assume there exists at least one Si G lCc{c), where 
Asi 7 ^ 0. 

Inductive Step Suppose we find a cycle c„+i G EC \ C 
where an edge exchange leads to Asi = 0. From 
Corollary 1^ subset C U c„+i will at least N + 1 unique 
sensors and the new cycle must introduce some sensor 
Sj ^ ICc- From Lemma 0 ( 2 ) As, 7 ^ 0. 


C. Proof of Lemma 

Proof: The matrix B G {—1,0,and a valid 
placement of of size /j.(G') — M — N + 1, therefore, Sat G 
{ —1, 0,-1-1}'^^^“^. Since the edges of G \ Ad maintain a 
spanning tree property, the graph has 1 connected component. 
Therefore the incidence matrix must be of rank of A — 1 czi. 
Finally, the matrix B'fj is a square matrix of size A — 1 and 
invertible. ■ 


3) Decoupling of Measurement along Cycle: We show that 
when an edge exchange occurs on a cycle, we need to only 
consider changes of flow values in JC{c). In the development 
/C(c) and edge exchanges are focused on cycles in EC. 

Proposition 2: A single edge exchange on c, with vertices 
{vq, ... ,Vm} results in a permutation of an uninterrupted path 
of the vertices. Therefore, if e 7 ^ e', pic, e) 7 ^ p{c, e'). 

Lemma 5: Consider a single edge exchange on c, e — e' 
where the following holds: 

PI Vs ^(c). As = 0, 

P2 Vs G E[c), As 7 ^ 0. 

Proof: (1) From Proposition]^ any sensor that measured 
a single vertex in p will measure all the vertices in p before 
and after the edge exchange. (2) From Proposition since we 
rearrange all the nodes yet keep a fixed edge to measure flows, 
all sensors in JC{c) will change values. ■ 

Consider the example in Figure | 8 (a)| where the cycle in 
consideration Ci will have a single edge exchange: for example 
Ae = (8,7) —>■ (1,2) where all Xi = 1. In this case, s = 
2 and and s' = —2. Since we assume that the magnitude 
and direction of the flow is measured, therefore there is no 
ambiguity. Any other edge will again lead to a new measured 
value. Since nodes {uo,...,uio} are always connected, all 
sensors outside of JC{c), for example a measurement on edge 
e = ( 2 , 10 ) will not detect an edge exchange. 

Remark 4: In condition (2) of Lemma we have As 7 I 0 
since we measure both magnitude and direction of flow. 


D. MAP Detector Structure 


Here we show how the general MAP detector rule can be 
evaluated for where we combine edge flows s, load forecasts 
X and candidate spanning tree T. 


T = argmaxPr(T | s, x) 

Tgt 

Pr (s, X I T) Pr (T) 

= argmax-— - 

Tgt Pr(s, x) 

= argmaxPr(s, x | T) Pr (T) 

Tgt 

= arg max Pr (s | x, T) Pr (x | T) Pr (T) 
Tgt 

= arg max Pr (s | x, T) Pr (x) Pr (T) 

Tgt 

= arg max Pr (s | x, T) 

Tgt 


(30) 

(31) 

(32) 

(33) 

(34) 

(35) 


Lines ( |30| ) - ( |3^ convert the MAP detector to a likelihood 
detector with prior weights. Line ([3^ conditions on the load 
forecast x. Since x does not depend on the outage hypothesis, 
(only s does), the term can be removed leading to ([35). 
Additionally, we assume a uniform prior over all hypotheses, 
however this does not have to be the case. 


E. Proof of Theorem 

Proof: We need to show that eq. d encodes the 
correct spanning tree: so f* = This can be done by 
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contradiction; Assume that the solution to f* = rep¬ 
resents some other spanning tree or even connected graph 
with some flow. This implies — i3]^s(72,x)) = 

— i3]^s(7i,x)). Since, exists, this im¬ 

plies that i?]^s(72,x) = i3]^s(7i,x). Which implies that 
s(7i,x) = s(72,x). However from Theorem 1, this is a 
contradiction. 

Finally, we must show that no subgraph of G' is uniquely 
distinguished from 7i. Such a subgraph occurs for any span¬ 
ning tree, where a removed edge of the co-tree is added to T. 
In such a case s(7i,x) ^ s(G",x) because: 

1) Repeat Lemma for not only edge exchanges, but the 
case where the cycle has no removed edge. In this case, 
we have again that [PI] and [P2] hold. 

2) The remaining proof can be repeated for not only all 
spanning trees, but subgraphs formed by these subgraphs 
formed by co-tree edge addition. 


F. Proof of Theorem 

Proof: We construct this test statistic 

Bh% - fr) = [fo(*i), • • •, fo(v)]^ (36) 

Where the matrix Bh S {0, where Bnik^ik) = 1 

for each edge in the co-tree of the particular spanning tree 
under hypothesis. 

The test statistic can be reduced to the following: 


Bh 


(37) 

Bh 

- {B]f~^)Blis{T, x)) 

(38) 

Bh 


(39) 


H 


Therefore, the new hypothesis test using x(T), is equivalent 
to the previous test using observed flow with some linear trans¬ 
formation H. The matrix is H is full rank in the \M.\ = f{G) 
case. Clearly if H is not square fA > H will not be full 

rank and the two tests will no longer be the same. ■ 



